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ABSTRACT 

A Chebyshev curve C{a, b, c, if) has a parametrization of the 
form x(t) = Ta{t)\ y{t) = Tb{t)\ z{t) = T^{t + ip), where 
a,b,c are integers, Tn{t) is the Chebyshev polynomial of 
degree n and (/? G R. When C(o, b, c, ip) has no double points, 
it defines a polynomial knot. We determine all possible knots 
when a, b and c are given. 
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1. INTRODUCTION 

It is known that every knot may be obtained from a poly- 
nomial embedding R ( [151 [5] ) . 

Chebyshev knots are polynomial analogue to Lissajous 
knots that have been studied by many authors (see [T] [51 
[8]|9][T3]). All knots are not Lissajous (for example the tre- 
foil and the figure-eight knot). In ID], it is proved that any 
knot K C R"^ is a Chebyshev knot, that is to say there exist 
positive integers a, b, c and a real such that K is isotopic 
to the curve 

C{a,b,c,ip) : x = Ta{t), y = Ti{t), z = T,{t + p), 

where T„ is the Chebyshev polynomial of degree n. This is 
our motivation for the study of curves C(a, b,c,ip), ip £ R. 

In |10) . the proof uses theorems on braids by Hoste, Zirbel 
and Lamm ([8] [13]), and a density argument (Kronecker 
theorem) . 

In [12], we developed an effective method to enumerate 
all the knots C{a,b,c,ip),Lp £ R where a = 3 or a = 4, a 
and b coprime. This method was based on continued frac- 
tion expansion theory in order to get the minimal b, on 
resultant computations in order to determine the critical 
values if for which C{a, b, c, ip) is singular, and on multi- 
precision interval arithmetic to determine the knot type of 
C{a,b,c,p). Our goal was to give an exhaustive list of the 
minimal parametrizations for the first 95 rational knots with 
less than 10 crossings. We obtained in [12] almost every min- 
imal parametrizations. For 6 of these knots, we know the 
minimal b and we know that the corresponding c must be 
> 300. 

In this paper, we develop a more efficient algorithm. It will 
give the parametrization of the 6 missing knots in [12] and 



allows to compute all diagrams corresponding to C{a, b, c, p), 
p (zH. One motivation is first to achieve the exhaustive list 
of certified minimal Chebyshev parametrizations for the first 
95 rational knots. Another is to provide a certified tool for 
the study of polynomial curves topology. 

Let us first recall some basic facts about knots. 

Knot diagrams 

The projection of the Chebyshev space curve C{a, b, c, p) 
on the (a;, t/)-plane is the Chebyshev curve C{a,b) : x = 
Ta{t), y — Tt{t). If a and 6 are coprime integers, the curve 
C{a,b,c,p) is singular if and only if it has some double 
points. It is convenient to consider the polynomials in Q[s, t, p] 

p _ T4t)-T4s) _ T4t + p)-T4s + p) 

t — S t — S 

We see that C{a, b, c, p) is a knot iff {(s, t), Pa{s, t) — Pb{s, t) = 
Qc(s, t, p) = 0} is empty. 

We shall study the diagram of the curve C{a,b,c,p), that 
is to say the plane projection C{a,b) onto the (x, 2/)-plane 
and the nature (under/over) of the crossings over the double 
points of C{a, b). There are two cases of crossing: the right 
twist and the left twist (see [IS] and Figure [T])). In [12| . 




Figure 1: The right twist and the left twist 

we showed that the nature of the crossing over the double 
point Ac, 13 corresponding to parameters {t — cos{a + P), s = 
cos(a — /?), a — /3 = ^), is given by the sign of 

D{s,t,p) = Qc{s,t,p)Pb^a{s,t,p). (2) 

D{s, t,p) > if and only if the crossing is a right twist. 

Note that the crossing points of the Chebyshev curve C{a, b) : 
X = Ta{t), y = Tb{t) lie on the (b-l) vertical lines T[,{x) = 
and on the (a — 1) horizontal lines Ta{y) = 0. We can repre- 
sent the knot diagram of C{a, b, c, r) by a billiard trajectory 
(see [10]). which is a pure combinatorial object. As an exam- 
ple, consider the knots 62 = C(4,5, 7,0), 52 = C(5,6, 7, 0), 
4i = C(3,5, 7, 0). We can represent their diagrams by the 
billiard trajectories in Figure [H 

When a = 3 or a = 4, we obtain the diagrams in the Conway 
normal form. In this case, the knot is rational and can 
be identified very easily by its Schubert fraction (see |15l 




52 52 4i 




Figure 2: Billiard trajectories 



111]). When fe > a > 5, the problem of classification is much 
more difficult. Nevertheless, the knowledge of the diagrams 
allows the computation of all the classical invariants, like 
the Conway, Alexander and Jones polynomials (see [15)). 

Summary 

Our goal is to compute all diagrams of C{a,b,c,(p), where 
a, b, c are given integers and <^ G R. From the algorithmic 
point of view, the description of the Chebyshev knots is 
strongly connected to the resolution of: 

Va,t,c = {Pa{s,t)=0, Pt{s,t)=0, Q,{s,t,<fi)=0}. (3) 

We first want to determine the set Za,b,c of critical val- 
ues for which the curve C{a,b, c, ifi) is singular. Because 
deg^ Qn = n — 1 and the leading term of Q„ is 2"~^rup"~^ , 
we showed in , that Va,b,c is 0-dimensional and has at 
most (a — 1)(6— l)(c— 1) points. We deduced that |.Ea6cl < 
i(a-l)(6-l)(c-l). 

Let Za.b.c = {vi, ■ • ■ , (^nj.The type of the knot C(a, b, c, is 
given by its diagram which is constant when tp is in {tpi, <^i+i), 
because the crossings do not change in this interval. In or- 
der to get all possible knots C(a, b, c, <p), we only need sample 
points ri in each {tpi,tpi+\) and to compute the diagram of 
C(a, 6, c, n). 

We can determine a polynomial Ra.b.c G Z[(/?] such that 
Za,b,c = Z{R). It can be defined by {R) = {Pa, Pb, Qc) fl QM 
and may be obtained with Grobner bases ([3]). In [T2j, we 
optimized the computation by an ad-hoc elimination based 
on the properties of the curves for a = 3 or a = 4. Grobner 
bases could fully be substituted by some resultant computa- 
tions in Z[X, ifi], the systems being generic enough. However, 
this leads to solve systems of very high degree. 

In the present paper we decompose the system by work- 
ing on some (real cyclotomic) extension fields. We show 
that the system ((3]) is equivalent to the resolution of ^ (a — 
!)(& — 1)[§J second-degree polynomials with coefficients in 
Q(cos ^, cos f , cos ^). This result is deduced from geomet- 
ric properties of the implicit Chebyshev curves. 

We show some properties of these extensions that allow to 
simplify the computations. We can represent the coefficients 
of the polynomials by intervals and certify the resolution. 
We then easily and independently obtain the roots of the 



second-order polynomials and the main difficulty becomes 
to compare them. A formal method would consist in com- 
puting their minimal polynomials over Q, which is equiva- 
lent to the resolution of (|3]). We use multi-precision interval 
arithmetic for coding the algebraic numbers cos as well 
as the solutions ip we get. If the two intervals are disjoint, 
the roots are distinct. If not, we can certify whether the re- 
sultant of the two second-order polynomials equals or not 
by Euclidean division. 

In section [21 we first describe the Chebyshev polynomials 
and the link between their factorizations and the minimal 
polynomials of cos ^ . This allows us to represent efficiently 
the elements of Q(cos cos ^, cos ^). Along the way, we 
give an explicit factorization of the Chebyshev polynomials. 

In section [3] we recall the definition of Lissajous curves and 
we give their implicit equations. We study the affine im- 
plicit Chebyshev curves Tn{x) = Tm(y) and show that they 
have [ t"^™-) J _|_ \ irreducible components, [ J being 
Lissajous curves. 

This allows us to deduce an explicit factorization of Ra.b.c as 
the product of second-degree polynomials Pa p -i, in section 

Bl 

We show in section [5l how to obtain Za,b,c, the set of roots 
of Ra.b.c, with their multiplicities. The general algorithm 
is described in [Bj This allows us to sample all Chebyshev 
knots C{a,b,c, </?), £ R, by choosing a rational number r 
in each component oi Yi — Za.b.c- 

In section [7] we find an exhaustive and complete list of the 
minimal parametrization for the first 95 rational knots. The 
worst case appears with the knot IO33 = C(4, 13, 856, 1/328), 
with deg -Ra.b.c = 15390. We discuss the efficiency of our 
algorithms and compare with those of [12[ . 

2. CHEBYSHEV POLYNOMIALS 

Chebyshev polynomials and their algebraic properties play 
a central role here. The curves we will study are defined by 
Chebyshev polynomials. The algebraic extensions we will 
consider are spanned by their roots and we need to know 
their factors. In this section we recall some classical proper- 
ties of Chebyshev polynomials. We will also show the link 
between their effective factorization in Q[t] and the minimal 
polynomial of cos ^ . 

The Chebyshev polynomials of the first kind are defined by 
the second-order linear recurrence 

To = 1, Ti =t, T„+i = 2tr„ -T„_i. (4) 

Til G 7j\t] and satisfies the identity Tn{cos9) — cos n9, and 
more generally r„ o Tm. = Tnm ■ We have 

n-l 
fc=0 

Let Vn be the Chebyshev polynomials of the second kmd 
defined by the second-order linear recurrence (the same as 
in (HI)) 

Vq = 0, = 1, K+i = 2tK - K-i. 



Vn G Z\t] and satisfies V„(cos6) — — . We have 

sm 

n-l 

V„ = 2'^-' Hit- COS 

and therefore Vd|Ki when d\n. Let us summarize some useful 
results in the following 

Lemma 1. We have the following properties: 

• Tl,{t) = T„(t) = ±1 

• T„{t) = ±1 ^ T^{t) = or t = ±1. 

• Tn{t) = y has n real solutions iff \y\ < 1. 

• Tn{t) = 1 has [^J real solutions. 

• Tn{t) = —1 has [-^^^J real solutions. 

Proof. From = nVn, we deduce that t i-> T„{t) is 
monotonic when \t\ > cos — , that Tn has n — 1 local extrema 
for tk = cos ^ and T„{tk) = (-1)'. □ 

Minimal Polynomial of cos — 

2i7r 

Let (n = e " ■ It is well known ([20j) that the degree of 
Q(Cn) is (p{n) where (p is the Euler function. Q(cos 221) = 
Q(("„)p|R and the minimal polynomial over Q of cos ^21 
has degree ^(p{n) when n > 1. Its roots are cos 2^ where 
k is coprime with n. Consequently, the minimal polynomial 
Mn of COS ^ has degree ^(p{2n), when n > 1. Its roots are 
tk ~ cos where (k,n) = 1, and k is odd. M„{—t) is the 
minimal polynomial of cos — . The leading coefficient of Af„ 

is 2'^(2n)/2^ 

Remark, cos G Q iff ^(p{2n) = 1 or n = 1, that is 
n = 1, 2, 3. In this case we get 2 cos ^ G Z. 

We deduce the following 

Proposition 2. LetPn be defined by Po = 1, Pi = 2i-l, 
P„+i = 2tP„ + Pn-i. Then we have (-1)"P„(-T2) = V2n+i 
and 

Pn^ U Ma (5) 

d|2n + l 

Proof. We have Po{-T2) = Vi, Pi{~T2) = -2r2 - 1 = 
— Vs. The sequences V2n+i and ( — l)"Pn(— Tb) satisfy the 
same recurrence formula: V2n+3 ~ 2T2V2n+i — V2n-i. Let 
d = 2m + 1 be a divisor of 2n + 1 and consider t = cos ^ = 
-cos 2^. We have (-l)"P„(t) = Vb„+i(cos = 0. 

Thus Md\Pn and we conclude using the fact that 

E degMd = I + l E <^(2d) =n = degP„. □ 

d\2n+l d\2n + l,d>l 

Lemma 3. We have Mjfc^ = Mm{T2h) tf m is odd. 

Proof. We have MmoT^k (cos ^) = and (2'=, m) = f 
so M2k^\M,n{T2k). We conclude since M^k^ and Mm(T2k) 
have same leading term. □ 

The relations between the minimal polynomial of cos ^ and 
the factorization of T^n^^_^—T^ri.^ is known ([20]). Formula 

([5]) together with Lemma [3] give also an algorithm to com- 
pute A/„. 

The number of factors of T„ is known ([7]). We give here 
the relation between the Chebyshev polynomials T„ and Vn 
and the polynomials M„. 



Proposition 4. Factorization of T„ and V„. 

We have the following factorizations in irreducible factors 

^2^(2^+1) = n (tlMd{T2,)y Md{t)Mdi-t) 

d\2m. + l \j=l / 

T2k(2m+1) = 2 n Md{T2k + l) 
d|2m+l 

where M„ is the minimal polynomial of cos ^ . 

Proof. The factorization of Vn is obtained by comparing 
its roots with those of A4d{H), when d\n. Let d be an odd 
divisor of n. We write n = 2^ ■ d\ ■ d, where di is odd. 
cos = cos 2k+id ^ of Tn SO M2k+id\Tn. We deduce 
the factorization by comparing the leading terms. □ 

3. CHEBYSHEV AND LISSAJOUS CURVES 

The following proposition will explain the notions of Lis- 
sajous and Chebyshev curves. 

Proposition 5. The parametric curve 

C : X = cos{at), y — cos{bt + ip), i G C, 

where a,b are coprime integers (a odd) and <^ G R admits 
the equation 

Tb(xf+Ta{yf -2cos{aip)Tb(x)Ta{y) - sm^{aip) = 0. (6) 

1. If a(p 7^ kiv, this equation is irreducible. C is called a 
Lissajous curve. Its real part is 1-1 parametrized for 
t G [0,27r]. 

2. If aip = fcyr, this equation is equivalent to Tt{x) = 
( — l)*ra(y). C is called a Chebyshev curve. It can be 
parametrized by x = Ta{t), y = { — l)''Tt,{t). 

Proof. Let {x,y) G C. Wehavert,(3;) = cos{bat), Ta(y) = 
cos{bat + aip). Let A = a(f, 9 = abt. We get Ta(y) = 
cos{e + X) so (l-cos^ 9) sin A = (cos 9 cos >^-Ta{y)f, that 
is {\ -Ta{x)f sin^ A = {Ta{x) cos \-Ta{y)f , and we deduce 
our Equation Q. 

Conversely, suppose that [x, y) satisfies Let x = cos(at) 
where t £ C. We also have x = cosa(t + 2|2l). We have 
Tbix) — cos 9. A — Ta{y) is solution of the second-degree 
equation 

— 2 cos{aip) cos 9 A — sin^ (oip) = 0. 

Consequently, we get Ta{y) = cos(^ ± aif) — Ta{cos(itbt + 
if)). We deduce that y = cos(±fet + + 2^), /i G Z. Chang- 
ing t by —t, we can suppose that 

X = cos at, y = cos{bt + ip + 2^). 

By choosing k such that kb + h = O(moda), we get x = 
cos at', y = cos{bt' + (f), where t' = t + 2^. 

If aip ^ 0(mod7r). Suppose that Equation ^ factors in 
P{x,y)Q{x,y). We can suppose, for analyticity reasons, 
that P{cos{at),cos{bt + (p)) = 0, for t G C. The curve C 
intersects the line y = in 2b distinct points so deg^ P > 2b. 
Similarly, deg^ P > 2a so that Q is a constant which proves 
that the equation is irreducible. 

licosaip = ( — 1)*^, the equation becomes ri,(a;) — (—l)*Ta(i/) — 
0. In this case the curve admits the announced parametriza- 
tion (see [51[TD] for more details). □ 



Remark. If a = & = 1, we obtain the Lissajous ellipses. 
They are the first curves studied by Lissajous (HI]). Let 
fi ^ (modTr). The curve 

£^ : x'^ + y'^ — 2cos{p)xy — sin^(/i) — 

is an ellipse inscribed in the square [—1,1]^. It admits the 
parametrization x = cost, y = cos(t + /i). 

The following notation will be useful. Let E^{x,y) = x^ + 
y^ — 2cos{ii)xy — sin^(/i) when fi ^ (modTr) and Eq = 
X — y, Ett = X + y. The Equation (|6]) is equivalent to 
Eaip{Tf,{x),Ta[y)) = 0. This shows that the real part of 
the curve C (Equation ((6|) is inscribed in the square [—1, 1]^. 
Using Proposition Owe recover the classical following result. 

Corollary 6. The Lissajous curve x — cos(at), y = 
cos{bt + ip), (aip ^ (modn) ) has 2ab — a — b singular points 
which are real double points. 

Proof, singular points of C satisfy Equation ((6)1 and the 

system 

Tl{x){Tb{x) - Ta{y)cosaip) = 0, 
T:,{y){Ta{y) - Tb{x) cos aifi) = 0. 

Suppose that Tl,{x) = Ti{y) = then T^{y) = T^{x) = 
1 and Equation |6] is not satisfied. Suppose that Tb{x) ~ 
Ta{y)cosaip = Ta{y) — Ti,{x)cosa(p — 0, then Ti,{x) = 
Ta{x) = and Equation |6] is not satisfied. We thus have 
either T^ix) = and Ta{y) ~Tt{x) cosaip — that gives (6 — 
1) X o real points because of the classical properties of Cheby- 
shev polynomials, or T^{y) ~ and Tt{x) — Ta{y) cosaip = 
that gives 6 x (a — 1) real double points. □ 

Remark. The study of the double points of Lissajous curves 
is classical (see [2] for their parameters values) . The study of 
the double points of Chebyshev curves is simpler (see I10|). 

Corollary 7. The affine implicit curve Tn{x) = Tm{y) 
has L^^J L'^^T^J + LfJ LtJ singular points that are real dou- 
ble points. 

Proof. The singular points satisfy either T„{x) — Tm{y) = 
1 or Tn{x) = Tm{y) = — 1 and we conclude using Lemma 

m □ 

Theorem 8. Factorization of Tn(x) — Tm(y). 

Let m = ad, n = bd, (a, b) — 1 and a odd. We have the 
factorization 

T„ix)-T^{y) = 2''-^l[Ck{x,y) 

fe=0 

where 

Cu{x,y) = E2akn {Tt{x),Ta{y)) 

d 

is the irreducible equation of the curve Ck '■ x — cos{at), y = 
cos{bt + 2^), given in Proposition\^ 

Proof. Let C be the curve T„{x) — Tm{y). We easily get 
Ck d C and Ck 7^ Cy . When fc = 0, Co admits the equation 
Tb{x)-Ta{x) = 0. If 2k = d, Ck admits the equation Tb{x) + 
Ta{y) — 0. In the other cases, the dominant term in x of 
Ck is 2^''"^J:^^ If d is even, we deduce that the dominant 

term of Yl Ck{x,y) is 2^'' ^'''a;^" and we get our result in 
this case. If d is odd, we get the same result. □ 



Corollary 9. Let d — gcd(a,6). The curve Th{x) = 
Ta{y) has +1 components. [^^J of them are Lissajous 
curves. 




Ti{x) = Tz(y) Tw{x) = Tz{y) 




Tio{x) = my) Tio{x) = Ts{y) 



Figure 3: Implicit Chebyshev curves 

Theorem |8] is particularly interesting when m — n = d and 
a = b = 1. In this case the curve Tn{x) = T^iy) is a union 
of ellipses and some lines. It will be useful for the determi- 
nation of the double points of Chebyshev space curves. We 
have 

Lfj 

= 2"-' n E2k. {s, t). (7) 
^ fe=i " 

The curve TuS^ll — '^"^'^-^ = has [^] irreducible compo- 
nents. Note that f 2fc7r and £21tv_ intersect at the point (t, s) = 

rt m 

(cos(-^ + ^),cos(^ — ^)) and its symmetric with respect 
to the lines s — —t and s = t. We recover the parametriza- 
tion of the double points of x = Ta{t) , y — Tb{t) that will be 
very useful for the description of Chebyshev space curves. 

Proposition 10 ([TOllIS]). Let a and b are nonnega- 
tive coprime integers, a being odd. Let the Chebyshev curve 
C be defined by x = Ta{t), y — Tb{t). The pairs {t,s) giving 
a crossing pomt are 

i = cos(^ + f), s = cos{^^^) 

where 1 < i < |(a - 1), 1 < j < & - 1. 

4. CRITICAL VALUES 

A polynomial Ra,b,c £ Z[(^] for which Za,b,c ~ Z{R) can 
be defined by {R) = (Pa, fs, Qc) H QM ^^'^ may be ob- 
tained with Grobner bases (0). 

Example. When a = 3, 6 = 4, c = 5, we find that 

Ra,b,c = (sotp-' + eo^" - 1) ■ 

(6400 ip^ - 3200 p'^ + 560 ip'^ - 80 (^^ -f l) . 
There are exactly 6 critical values that are symmetrical 




t - s 




Figure 4: Double points in the parameters space 
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1/3 = .117, .128, .590 
Figure 5: Pa = 0, P4 = 0, Qs = 

about the origin. For these values of tp, the curve Q5 (s, t, ip) = 
0, which is translated from the curve Psis, t) = by the vec- 
tor {tfi, if), meets the points {P3 = 0, P4 = 0}. 

In this part, we use use the properties of Chebyshev curves 
obtained in section [S] We give an explicit formula for the 
polynomial Ra,b,c as a product of of univariate polynomials 
of degree 1 or 2 with coefficients in Q(cos ^, cos f, cos ^). 

Proposition 11. Let a,b be nonnegative coprime inte- 
gers and c be an integer. Suppose that a is odd. Let Ra^b,c{f) 
be the polynomial 

a~l 



2 6-1 



n llQcicOs{i; + j;)T:,COsi 



(8) 



Ra,b,c £ and C{a,b, c, If) is singular iff Ra,b,c{^) ~0. 

Proof, ip G Za,b,c iff there exists (s, t) such that Pa{s, t) = 
Pb{s,t) = and Qc{s,t,Lp) — 0. This conditions are equiv- 
alent to have t = cos(^ -|- ^) and s = cos('2^ - ^) and 
Qc{s,t,(p) — 0, for some 1 < i < ^^-ji and 1 < j < & — 1, 
from Proposition 1101 

Qc{s,t,if) is a symmetric polynomial of Z[iy9][t,s]. Let — 
cos(Qi + Pj), t = cos(a. 



s + t = 2cosQiCos/3j and st = cos^ 4- cos^ /3j — 1, we 
deduce that Qc{s,t,ip) belongs to Z[</3, cos Qfi][cos 

6-1 

Ri = ]^Qc(cos(ai -|-/3j),cos(Qi 
3=1 

belongs to Z[(/3,cosQi] because the cos /3j are the roots of 
Vb G Z[t]. From Qc{ — s,—t,ip) — Qc(s,t,~(p) we deduce 

a-l 

2 a-l 

that n Ri{-ip)Ri{ip) = n Riiv) G Z[</j]. Wc thus have 



R: 



a,b,c 



G TiVp] and so it is for Ra. 



b,c. 



□ 



Let s = cos(a + /?), t = cos(a — /3). Using Theorem |B] and 
Formula ((T)), we get 



Qc(s,t,ip) ^2''-^Y[E2kj^{s,t). 



k-l 



Let us consider Pa s — — 5 — E2~,{s ip,t-\- ip). For 7 7^ 
' ' sin 7 

-Pq,/3,7 is 

2 „ ^ (cos^ Q — cos'^ 7)(cos^ /3 — cos^ 7) 
(p^ + 2vp COS a COS ^ -f -5^ '^^ ^ 



2 ' 



sin^ 7 



and 



-P n 2L = + COS a COS /3. 



We therefore obtain 



Qc(cos(a -I- /3),cos(a - /3), (p) = t^- (vs) 



with if = 2"-^ n 2sin ^ = c2'=~\ We get therefore 

c I g — 1 

1 LjJ 2 iJ-1 

Ra,bA'p) = if 2 (-DC'-i) n n n 2i M (v')- 

fc = l ,=1 j = l a ' i, ' c 

We have written Ra,b,c as the product of second or first- 
degree polynomials Pa,/!,-, in Q(cos cos ^, cos ^). 

5. COMPUTING THE CRITICAL VALUES 

Our strategy consists in first computing separately the real 
roots of each Pa,p,~t and then combining these roots to get 
those of Ra,b,c- A straightforward approach would be to use 
interval arithmetic to approximate the various trigonometric 
expressions, but this would fail when Ra,b,c has multiple 
roots, unless we cannot ensure if some discriminant or some 
resultant are null. 

5.1 Real roots of Pa ,3,7 

Let Q = ^, /? = ^ and 7 = with 1 < i < 1 < 

a'*^ b ' c — — 2' — 

j<b-i, i<k< L^j. 

If 7 = -J , the unique root of P^ ^ ir is — cos a cos /3. If 7 7^ § , 

11... r 7-. -A 9/1 sin^ asm^ B 
the discnmmant of Pa,p,-y is 4 cos 7 1 1 — 



has the same sign as 



sm 7 



It 



^, ft- = ^ and s 



Pj). From 



sin^ 7 — sin^ a sin^ /? (9) 

The knowledge of the sign of ((9)| then gives explicit formulas 
for the real roots of Pa,p.-y. 



5.2 Multiplicity of 

Proposition 12. The multipUcity of ip = in Ra,b,c is 

V((6,c)-1)+L|j((a,c)-1). 

Proof. We have to examine whenever (^3 = is a root of 
Pa, 13,1 where a = /3 = ^ and 7 = Here a is odd so 
cos a 7^ 0. Thus, = is a root of ^ ir if and only if 

cos a cos /3 (10) 

is null and when 77^^, (^ = is a root of Pa,/3,-y if and only 
if the following expression is null 

(cos^ a — cos^ 7)(cos^ /3 — cos^ 7). (11) 

• If7 = /3=-|,</7 = 0isa root for i = 1, . . . , 2-^. 

• If77^-|,(y3 = 0isa root of -Pa, ^,7 if and only if 
sin^ 7 = sin^ a or sin^ 7 = sin^ /3, that is ic = ka or 
jc — kb or (b — j)c — kb. The root (p — is obtained 
for i = X-r^, k = X-r^, X = l,...,-^^4^ and it 

(a,c) ' (ci.c) ' ' ' 2 

is double when P — ^. It is also obtained for j — 
fJ-{ET)' ^ = ^'I^' M = l.---,(6,c) - 1. We obtain 
L|jaa,c)-l) + ((6,c)-l)(a-l)/2. 

We thus obtain the result. □ 

Remark. We find that is not a critical value if and only 
if a, b and c are pairwise coprime integers. This result was 
first proved by Comstock ([3], 1897), who found the number 
of crossing points of the harmonic curve parametrized by 
x^Ta{t),y = Tt{t),z^T,{t). 

5.3 Non null multiple roots of Ra,b,c 

It may happen that Ra.b,c has multiple root if. Several 
cases may occur. 

► Pa,i3,-y has a double root if and only if Disc {Pci,i3,-t) = 0, 
that is to say sin^ 7 — sin^asin^/3. The double root is 
ip — — cos a cos p. 

► Pa,i3,ji and Pa, ^,72 have a common root. In this case 
Pa,i3,-fi = Pa, /3, 72 1 that is to Say 

(sin^ 71 — sin^ 72)(sin^ 71 sin^ 72 — sin^ a sin^ /J) (12) 

is null. 

► Pai,/3i,7i and Pci2,P2,i2 have a common root. 
The first two cases are related to the equation 

sin riTT sin r27r = sin ravr sin r47r (13) 

where € Q. All the solutions of Equation (|13p are known 
(see [IS]). There is a one-parameter infinite family of solu- 
tions corresponding to 

sin ^ sin6' = sin | sin(-| — |), 

and a finite number of solutions listed in [16]. We deduce 
from a careful study of the Equation (|13p : 

Proposition 13. Let a ^ —, = ^ and 7 = — , where 

a ' i^ b ' c ' 

(a, b) = 1 and a is odd. Pa, ,9, 7 has a double root iff P — ^ 
and sin 7 — sin a. In this case, the double root is if — 0. 

and 



Proposition 14. Leta=if,p^if and 71 = 
72 = -^^j where {a,b) — 1 and a is odd. Then Pa, ^,71 and 
Pa, 13,12 have a common root if iff there are equal and 

1. sin Of = sin 71, sin/3 = sin 72. 

In this case, the roots are if = Q and if = —2 cos a cos /3. 

2. sin/3 = i, sin 71 = sin ^a, sin 72 = cos ^a. 

In this case the common roots are ip = — cos(a ± -g-)- 

3. sin 72 = sin 71 . 

► In case when Oi / 0,2 or /3i / /32, Pai„Si,7i and Pa2,P2,ii 
have a common root if Res ^(Pai,/3i,7i , Pa2,/32,72) ~ 0- This 
resultant can be expanded and its sign is the one of: 

((cos^ Qi — cos^ 7i)(cos^ /3i — cos^ 71) sin^ 72 — 

(cos^ a2 — cos^ 72)(cos^ /32 — cos^ 72) sin^ 71)^ 
— 4(cos Qfi cos P\ — cos a.2 cos Pi) sin^ 71 sin'^ 72 x 
((cos^ Qi — cos^ 7i)(cos^ P\ — cos^ 71) cos a2 cos Pi sin^ 72 — 
(cos^ Q2 — cof? 72)(cos^ Pi — cos^ 72) COS cvi COS /3i sin^ 71) . 

(14) 

It would be interesting to get an arithmetic condition as- 
serting that this resultant is null. 

5.4 Computing the diagrams 

Let If £ H. ip may be a rational number r G Q — Za,b,c 
or an algebraic number given by a polynomial whom it is a 
root and an isolating interval. The main step is the com- 
putation of the crossing nature at the double point Aa,p 
corresponding to parameters (t = cos a + P,s = cos a — /3), 
where a = — , P — There are two cases to consider. 

a ' ^ b 

1. We know the roots ipi < ■ ■ ■ < ipm of Qc{s, t, ip). 

\i f < ipi then n = Q otherwise let n = maxjfc, if > 
ipk}. We have sign [Qc{s,t,ip)) = (-1)". 

2. We do not know the roots of Qc{s, t, ip). 

Wc compute Qc{s,t,ip) using the recurrence formula: 

Qo = 0, gi = 1, ga = 25 + 4v5, 

Q3 = ^4:T+I2if S + A + 12 if'^ -3. 

g„ + 4 = 2{S + 2lf){Q,^ + 3 + Qn + l) 

-2{2ip'^ + 2T + 2ipS+l) Q„+i - g„. 

where S = s + t = 2 cos a cos /3 and T = st — cos^ a + 
cos^/3-1 (see [12]). We work formally in Q [it, •!;]/{ Af, iV) 
where M, N are the minimal polynomials of u — cos a, 
V = cos p. 

The sign of the crossing is 

D{s,t,lf) = Qc(s,t,ip)Pb-ais,t,ip) 

= (_l)'+^sini^sinifig,(s,t,^) 
= (-l)'+^+L-J + L-JQ^(,,t,^). 

6. THE ALGORITHM 

We want to compute all the real roots pi < . . . < ipn of 
Ra,b,c that factors in ^(a— 1)(&— 1) polynomials Pa^.0j,j^. 
We precisely want non overlapping intervals [am,bm] for 
these roots in order to chose sample rational points ro < ai, 
hi <ri < Oi+i, b„ < r„. 

At some stages, one may need to compute the sign of 
Disc (Pa,;3,7) (expression (j^J) or Res (Pai,/3i,7i , Pai,;3i,7i) (ex- 
pressions (|12|l and (|14|l ') in order to decide whether two 
roots are distinct or not. This information is required for 



two reasons. We first want to be sure that we get all the 
roots and secondly, we will need to know all the roots of 
Qc(cos(Q:i + Pj), cos (a,; — Pj), (fi) with their multiplicities in 
order to determine the nature of the crossing over the cor- 
responding double point in the diagram (section 15.41 

The signs of (|10|) and may be evaluated by simple 

arithmetic considerations on a, /3, 7. 

Isolate and Refine. A very first step is to get accurate 
isolating intervals with rational bounds for cos ai , cos /?-, and 
cos7fe to perform interval arithmetic for the real roots of 

Such intervals can be computed by performing algorithms 
based on Descarte rule of signs (see for example 18^) on the 
used Chebyshev polynomials V„. Algorithms like in [IS] can 
easily solve such polynomials for very high degrees (several 
thousands) with a large accuracy. The computation of the 
required isolating intervals can then be performed as a pre- 
processing for the global algorithm. 

From now, we denote by Isolate(P, acc) the function that 
isolates the roots of a univariate polynomial P with rational 
coefficients by means of intervals with rational bounds for a 
given accuracy acc (maximal length of the intervals). This 
function provides non overlapping intervals that contains a 
unique real root of P (and such that each real root of P is 
contained in one of the intervals). 

Note that if more accuracy is required for some intervals, 
it is easy to refine them from the isolating intervals provided 
by the function Isolate: just evaluating P at some points, 
without running again the Isolate function with a higher 
value for acc). We name by Refine(7,P, acc) the function 
that decreases the length of the interval / to get an accuracy 
< acc, knowing that the interval isolates a real root of P. 

IsolateP. Thanks to Proposition 1131 one can compute the 
roots of Pa,i3,-t with an appropriate accuracy, using multi- 
precision interval arithmetic for the evaluations. We will 
use the function IsolateP(Q:, /?, 7, acc) that returns a (pos- 
sibly empty) list of (a, /?, 7, [u, v]) corresponding to isolating 
intervals [u, v] for the roots. 

SignTest. When two isolating intervals [mi, vi] and [u2,«2] 
corresponding to (ai,/?i,7i) and ((a2,/32,72) are such that 
wi] ni'^a, V2] 7^ 0, we first use a filter (named SignTest in 
the sequel) which consists in using multi-precision interval 
arithmetic for the evaluation of Res <p(Pai,/3i,-yi , Pa2.fe.72 ) 
(expressions 1 12|) and (|14|) . 

Thanks to Proposition 1141 we know by arithmetic consid- 
erations when (|12|) is null. We know also the corresponding 
common roots and we change [11^,11;] to [mi,wi] n[w2,W2]. 

Expression (|14|) is 

p = {id - cDid - c|)(i - d) - 

(C|-C|)(C|-C|)(1-C|))' 
-4(CiC3 - C2C4)(1 - C|)(l - Ci)x 

{{c! - cl){cl - ci)C2C4{i - cl) - 

(pl - Cl)(Cl - Cl)CrCz(\ - Ci)) , 
where Ci = cosQi, C2 = cos 02, C3 — cos/3i, C4 — cos/32, 
C5 = cos 71, Ce = cos 72. 

Given isolating intervals with rational bounds that con- 
tain the values of the required d, i — 1, . . . , 6, the func- 
tion SignTest (ai , 02 , /3i , /32 , 71 , 72) straightforwardly evalu- 
ates P. If the resulting interval is [0, 0] or do not contains 0, 
one can decide the sign of the input, otherwise, the function 
returns FAIL. 

FormalNullTest. In case of failure of SignTest, one has 



to decide if the input is null or not, which is the goal of the 
function FormalNullTest we now describe. 

Let us write a, = ^,02 = -^,^2^ ^,7i = 

1 72 ~ ■ Let m be the smallest common multiple of 
ai , a2 , foi , 62 , Cl and C2 . According to the definitions of T„ , we 
have Ci — T„.(cos ^). Since Mm is the minimal polynomial 
of cos — , the expression P(Ci, . . . , Ca) is null if and only if 
P(r„,,'?. . ,r„j = in Q[t]/{Mmit)). 

DoubleTest. Our function first performs the SignTest. If 
it returns an interval with bounds of same sign, then the 
sign of the tested expression is the sign of the two bounds 
of the interval. Otherwise, we run the FormalNullTest. If 
this test returns then the expression is null. Otherwise, we 
decrease the lengths of the intervals that represent the values 
of cos ^ by calling the function Refine until the SignTest 
does not FAIL (the fact that the sign of the expression to 
be tested is known not to be ensures that this process will 
end). 

The global algorithm. We proceed in three steps : 

(0) We isolate the roots of some Chebyshev polynomials 
using the Isolate black-box with an arbitrary accuracy. 

(1) We compute separately the roots of the Pa,i3,-y by using 
IsolateP. 

(2) We then consider the list of these roots and observe 
carefully the overlapping intervals. For any pair of overlap- 
ping interval, we decide whether corresponding resultants 
are null or not using DoubleTest. If the corresponding roots 
are equal then we change their isolating intervals by taking 
their intersection. 

From these disjoint intervals with rational bounds, we straight- 
forwardly get the roots with their multiplicities. We thus 
deduce the sample points ro, . . . ,r„ we need. Furthermore, 
for each ai = ^1 Pj ~ '^1 we know the roots with their 
multiplicities of Qc{t, s, Lp), where t = cos(ai + Pj) and 
s = cos(ai — Pj). This information is helpful for knowing 
the crossing nature at the point Aa-^p^ (section l5.4p . 

7. EXPERIMENTS 

In the appendix of [12] , we gave parametrizations of every 
rational knot as C(3, b, c, (p) and C(4, fe, c, (p) where (&, c) were 
minimal for the lexicographic order (c < 300). For 6 knots 
we knew the minimal b and that c > 300. With the method 
we developed here, we recover all the minimal parametriza- 
tions we gave in [T^ but also for the 6 missing knots. The 
following knots admit the parametrizations: 

95 = C(3, 13, 326, 1/85), IO3 = C(4, 13, 348, 1/138), 
IO30 = C(4, 13, 306, 1/738), IO33 = C(4, 13, 856, 1/328), 
IO36 = C(3, 14, 385, 1/146), IO39 = C(3, 14, 373, 1/182). 

For example, one deduces that there is no parametrization 
of 95 as Chebyshev knots with (a, 6, c) <iex (3, 13, 326). 

P3,i4,385 has degree 4992. It has 2883 real roots. All are 
simple roots except that is of multiplicity 6. 

P4,i3,856 has degree 15390 and 9246 real roots (0 has mul- 
tiplicity 18). We get 2050 non trivial knots, 83 of them 
are distinct, and 63 have less than 10 crossings. The to- 
tal running time — critical values with their multiplicities, 
sampling of 1442 values, computing knot invariant — was 
450" (Maple 13, on Laptop, 3Gb of RAM, 3MHz). 

Outside the intrinsic combinatorial aspects of the problem, 
the complexity of our algorithm essentially depends on the 



FormalNullTest. In the worst case d = abc and degM^ = 
|(a — l)(c— 1), when a, b, c axe prime integers, the most 
difficult computation consists in deciding if the expression 
I14l is null or not which is equivalent to testing if a univariate 
polynomial of degree at most 4d is null modulo Md or not. 

These computations can be speed up a lot since they can 
be performed modulo a prime integer: all the considered 
polynomials have a power of two as leading coefficient and 
we just need to test if one polynomial is null modulo another 
one. 

In this challenging experiments, we never had to run the 
FormalNullTest, the SignTest being always sufficient, thanks 
to the filters given by propositions 1131 and 1141 and to a good 
(experimental) choice initial choice of accuracy when com- 
puting the prerequisites running the Isolate algorithm. 

8. CONCLUSION 

The method we developed in this paper allows us to com- 
pute Chebyshev knot diagrams for high values of a, b and c. 
Our experience with small a and b shows that the difficult 
cases (multiple roots of Ra,b,c we found) were predictable. 
There are certainly some specific reasons connected with 
arithmetic properties and the structure of cyclic extensions. 

The main difference with the algorithm described in [12] 
and the computation of Ra,b,c as a polynomial of degree 
^{a — l){b — l)(c — 1), is that it came as a resultant of a 
polynomial of degree (c — 1) in {X, tp) and a polynomial of 
degree f (a — l){b — 1) in X with coefficients in a unique 
field extension. The example described in this section can 
be considered as the extremal case, in terms of degree, to be 
solved using methods from the state of the art when running 
[12j while it can be solved in few minutes with the method 
proposed in this article. 

From the point of view of knot theory, it is proved in [llj . 
that rational knots with A'^ crossings can be parametrized by 
polynomials of degrees (3, b, c) where b + c < 3N which is far 
better that the results we obtain here. But we challenged to 
give, as it was done with Lissajous knots in [1], an exhaustive 
and certified list of minimal parametrizations. We consider 
that it might be one step in the computing of polynomial 
curves topology. 
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